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ABSTRACT 

Motivated by suggestions of 'cosmic downsizing', in which the dominant contribution to 
the cosmic star formation rate density (SFRD) proceeds from higher to lower mass galax- 
ies with increasing cosmic time, we describe the design and implementation of the Redshift 
One LDSS3 Emission line Survey (ROLES). This survey is designed to probe low mass, z^l 
galaxies directly for the first time with spectroscopy. ROLES is a X-selected (22.5 < A'ab < 
24.0) survey for dwarf galaxies [8.5 ^log(M*/M )^ 9.5] at 0.89 < z < 1.15 drawn from 
two extremely deep fields (GOODS-S and MS1054-FIRES). 

Using the [O II]A3727 emission line, we obtain redshifts and star-formation rates (SFRs) 
for star-forming galaxies down to a limit of ~ O.^M^yr . We present the [O II] luminosity 
function measured in ROLES and find a faint end slope of a. faint ~ —1-5, similar to that 
measured at z^O.l in the SDSS. By combining ROLES with higher mass surveys (GDDS 
and ESO GOOD-S public spectroscopy) we measure the SFRD as a function of stellar mass 
using [O II] (with and without various empirical corrections), and using SED-fitting to obtain 
the SFR from the rest-frame UV luminosity for galaxies with spectroscopic redshifts. Our 
best estimate of the corrected [O II] -SFRD and UV SFRD both independently show that the 
SFRD evolves equally for galaxies of all masses between z~ 1 and z^0. 1 . The exact evolution 
in normalisation depends on the indicator used, with the [Oil] -based estimate showing a 
change of a factor of w2.6 and the UV-based a factor of «6. We discuss possible reasons for 
the discrepancy in normalisation between the indicators, but note that the magnitude of this 
uncertainty is comparable to the discrepancy between indicators seen in other z~l works. 
Our result that the shape of the SFRD as a function of stellar mass (and hence the mass range 
of galaxies dominating the SFRD) does not evolve between z~l and z^O.l is robust to the 
choice of indicator. 
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1 INTRODUCTION 

Over the last decade or so, dozens of studies have measured the star 



formation rates of statistical samples of galaxies (e.g. , Lillv et al 



19961 : lMadauetal. |[l996; Steide l et alj|l999t iHopkins & Beacon: 
2006; Redd v & Steide| | 2009l). recently out to redshifts as high as 
z~6 and beyond (e.g. jBouwens et alj |2007). However, the major- 
ity of these are limited to the most actively star-forming and/or 
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the most massive galaxies. In order to construct a complete cen- 
sus of star-formation, it is necessary to probe to low enough star- 
formation rates (SFRs) to ensure that the star formation rate density 
(SFRD) has converged. Star-forming galaxies may exhibit a rela- 
tively narr ow sequence in st ar-formation rate as a function of stel- 
lar mass dNoeske et alj|2007r) . If this is the case, then the bulk of 
the high SFRs will be confined to the highest stellar mass objects. 
However, since low mass galaxies are much more numerous than 
their higher mass counterparts, it is an open question in what mass 
objects the majority of the SFRD resides. Even at z~l, when the 
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Universe was 40% of its current age, only the most massive galax- 
ies (los(Af»/Ma) >10.5) are routinely explored by spectro scopic 
observations (e.g. JCowie & Bareer l2008l:lDamen et al.ll2009h . The 
Gemini Deep Deep Survey (GDDS. I Abraham et alj|2004l) recently 
managed to obtain spectroscopy of much lower stellar mass objects 
(log(Af* /Mq ) > 10) but this required large amounts (~30 hours) of 
integration on an 8-m telescope. Lower masses can be investigated 
using photometric redshifts, but these either lead to larger uncer- 
tainties on the SFR for each using SED-fitting, si nce the redshift, 
SFR and dust must be fitted simultaneously (e.g.. lMobasher et"al] 
2009); or require statistical stacking of man y objects using an- 
other SFR indicator such as radio continuum jDunne et al. 2009; 
IPannella et al]|2 009). The drawback with the stacking technique is 
that systematic errors in the stacked sample may bias the results, 
and only the properties of an average galaxy can be studied, not the 
variation between galaxies. 

The combination of SFR and stellar mass data from sur- 
veys conducted at multiple epochs will provide important infor- 
mation on the evolution of star-formation and the build-up of stel- 



lar mass. For example , many recent studie s (fNoeske et al.l 



Cowie & Barger |2008l ; iDunne et alj 120091 ; IPannella et al. 



2007; 



2009; 



Damen et al . 2009, an d others) are finding evidence for the 'cosmic 
downsizing' of lCowie et al.l dl996h in which the dominant contribu- 
tion to the SFRD proceeds to progressively less massive galaxies as 
the Universe ages. However, such results have been based on direct 
observations of only the most massive galaxies at higher redshifts. 
Thus, it is an open question whether low mass galaxies have always 
had similar SFRs at all times (and thus 'downsizing' would simply 
be described by the decrease in the SFRD of high mass galaxies 
with cosmic time), or whether low mass galaxies show some dif- 
ferent evolutionary behaviour. In the latter case, it is conceivable 
that, at some point in the past, low mass galaxies may have dom- 
inated the universal SFRD. Thus, a useful way to directly explore 
these scenarios is by looking at the SFRD as a function of stellar 
mass, down to low stellar masses (log(M*/M0)~9), and how this 
evolves with redshift. 

Ideally, one would like to obtain a stellar mass-selected survey, 
complete to low SFRs and repeat this at multiple epochs/redshifts. 
To this end, we have designed the Redshift One LDSS3 Emis- 
sion line Survey (ROLES), using guaranteed time o n the LDSS3 
spectr ograph as part of the LDSS3 instrument project. Ipavies et ail 
(2009) (hereafter Paper I) used this survey of low mass galaxies 
at z~l (a subsample of dat a from the current paper) plus higher 
mass galaxies from GDDS ( jjuneau et alj|200 5l . hereafter J05) and 
compared this with a local measurement from the Sloan Digital 
Sky Survey (SDSS). In Paper 1, we found a suggestion that the 
[O II] luminosity density/SFRD at z-~l began to turnover below a 
mass of around log(M,/M )~9.5, and thus the bulk of the SFRD 
is in higher mass objects. Paper 1 presented initial results from a 
subsample of three masks in one of our pointings. Here we present 
all the data from the same redshift bin for both our fields (three 
different pointing positions), resulting in five times the number of 
z~l galaxies due to the greater areal coverage and higher spectro- 
scopic completeness . We also employ a new empirical correction 
dGilbank et all20ld hereafter G10) to accurately correct the mass- 
dependence of [O II]-SFRs; SFRs from SED-fitting; new higher 
mass samples from other spectroscopic surveys; and an improved 
local comparison sample built from the SDSS. 

This paper presents the survey design and method and re- 
sults using the complete dataset from the 0.89 < z < 1.15 
component of the survey. The layout of this paper is as follows. 
§2 describes the survey design, observations, and data reduction. 



§3 describes the analysis of the catalogues so produced and com- 
parison samples from other z~l surveys. §4 shows and discusses 
the results from the [O II] luminosity function and star formation 
rate density, and §5 concludes. All magnitudes are quoted on the 
AB system unless otherwise stated, and we assume a cosmology 
(h, Q,m, ^a) = (0. 7, 0.3, 0.7). Throughout, we convert all quanti- 
ties to those using a lBaldrv & Glazebrookl fe003h (hereafter BG03) 
universal initial mass function (IMF). 



2 METHOD 

2.1 Survey Design 

ROLES (the Redshift One LDSS3 Emission line Survey) is de- 
signed to probe the SFRD in galaxies with stellar masses much 
lower than previously studied at redshifts of order unity. In order to 
do this efficiently, we adopt a novel survey strategy. We utilise fields 
with very deep Tf-band photometry and photometric redshifts in 
order to pre-select galaxies which are most likely low-stellar mass 
systems at z~l. These galaxies are then followed up with multi- 
object optical spectroscopy in order to obtain both the redshift and 
star-formation rate. To increase efficiency even further, we use a set 
of custom band-limiting filters to restrict the wavelength range (and 
hence the corresponding redshift range for a given rest-frame wave- 
length) of our spectra and provide us with a high sampling density. 
We select A"-fahrQ (22.5 < K < 24.0) targets. Since securing ab- 
sorption line redshifts for galaxi es this faint is extrem ely challeng- 
ing (c.f. the brighter sample of lAbraham et alj|2004l . GDDS) and 
would require prohibitive amounts of telescope time, we choose to 
use relatively modest exposure times (~ 4 hours on a 6.5-m tele- 
scope) and specifically target the [O II] emission line within our 
redshift window. With this approach, we will not obtain redshifts 
for galaxies without emission lines (to our flux limit), but these do 
not contribute significantly to the SFRD. In this way, we obtain a 
mass-selected sample, complete to a given SFR limit (to the extent 
that [O II] traces the SFR). 

Fields with both the depth and wealth of multi-colour imag- 
ing we require for this survey are currently scarce. For ROLES, 
we use photometric catalogues and photometric redshifts from 
two imaging surveys in two different fields to select targets for 
follow-up spectroscopy: the Great Observatories Origins Deep Sur- 
vey (GOODS) region o f the Chandra Deep Field South (CDFS, 
e.g., IWuvts et alj |2008|). and the Faint Infra-Red Extragalactic 



Survey (FIRES, e.g.. iForster Schreiber et alj|2006l) . We note that 



FIRES actually comprises two different fields: the Hubble Deep 
Field South (HDF-S) and a region containing the z=0.83 cluster 
MS 1054-03. In ROLES we have only followed up the MS 1054-03 
area (due to its significantly wider area than the HDF-S), but here- 
after we refer to the MS 1054-03 field as "FIRES". The layout of 
our fields is shown in Fig.Q] 



2.2 Sample selection & Spectroscopic Observations 

For our LDSS3 spectroscopy, we obtained a set of four custom- 
band limiting filters to observe galaxies in well-defined redshift 
bins. In this work, we only consider the KG750 filter which cov- 
ers approximately (750 ± 50) nm. Observations using the KG650 
[(650±50) nm] are ongoing and will be presented in future work. 



1 all our magnitudes are on the AB system unless otherwise noted. K^b 
-f^Voga + 1-87 
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Figure 1. Layout of the ROLES fields: CDFS (left) and FIRES (right) on the same scale. CDFS comprises two pointing positions centred on 03 32 m 27.6 s , 
-27 d 45 m 00 s and03 h 32 m 28.8 s , -27 d 52 m 12 s of roughly 8 arcmin diameter, set by the field of view of LDSS3. FIRES is a single pointing of 5.5x5.3 arcmin, 
set by the limit of the deep photometry available. Points indicate galaxies with 22. 5< K < 24.0 and filled circles show our detections of z~l emission line 
galaxies. 
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Figure 2. The transmission curve of the KG750 filter. The dashed lines 
at the half-maxima denote the minimum and maximum wavelengths used, 
7040A-8010A. 



The transmission curve of the KG750 filter is shown in Fig.H Half- 
maxima occur at 7040 A and 8010A and these are adopted as the 
wavelength limits. In targeting the [O II] line at 3727 A, these lim- 
its correspond to a redshift range of 0.889 < z 1.149. 

LDSS3 can perform multi-object spectroscopy (MPS) in nod - 
and-shuffie (N&S) mode dGlazebrook & Bland-Hawthornl l200ll) . 
which allows extremely accurate sky-subtraction thus allowing 
lower flux limits to be reached than would otherwise be achiev- 
able with traditional MOS. This is particularly important in the 



wavelength window of interest here, which is populated by nu- 
merous bright sky lines. The LDSS3 medium red grating (300 
lines/mm was used, giving an average dispersion of ~ 2.7A/pixel) 
and 0.8" wide slits which, with the plate scale of 0.189"/pixel, 
gives a resolution of 1 1 .4A FWHM. The N&S observations utilised 
3"long slits, observing each target galaxy for 60s at one end of the 
slit and then nodding the telescope 1.2"and observing for another 
60s at a second position within the slit. These two positions are re- 
ferred to as A and B for convenience. By nodding within the slit, 
the galaxy is observed for the total time the detector shutter is open. 
Charge is shuffled along the detector between the A and B observa- 
tions by 16 pixels (corresponding to the slit width of 3"), resulting 
in two exposures for each galaxy, stored in different locations on 
the CCD. Because the telescope was nodded between exposures, 
the B exposure contains observations of the sky at the location of 
the galaxy in the A exposure and vice- versa. Thus sky-subtraction 
can be achieved simply by subtracting the A observation from the 
B. This process is described in more detail in the next section. 

With the above parameters, we can place an average of al- 
most 200 objects in each mask over the ~8.2 arcmin diameter 
circular field of view of LDSS3. To select targets, galaxies were 
prioritised by assigning a weight to them based on their A'-band 
magnitude, photometric redshift and its error. Galaxies in the range 
22.5 < K ^ 24.0 were primarily targeted. Photometric redshifts 
and confidence intervals kindly provided by B. Mobashe r and T. 
Dahlen (CDFS) and from iForster Schreiber et alj d2006L FIRES) 
were used to give higher weighting to those galaxies whose error 
bars placed them within (or overlapping with) the redshift window 
0.889 < z ^ 1.149. Once weights were assigned to each galaxy, 
slits were randomly allocated to objects, sorted by priority. Due to 
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Table 1. Summary of the observations. The two pointing positions for 
CDFS are denoted CDFS.l and CDFS.2 (northern and southern pointing, 
respectively). Total exposure times include only data of acceptable quality 
used in the analysis. Results based on masks indicated with "*" were pre- 
sented in paper 1 , but these have been re-reduced with the present method 
for uniformity. 



Mask 


Total exposure 
time Iks 


Field 


maskl 


31.2 


FIRES 


mask2 


24.0 


FIRES 


mask4* 


21.6 


CDFS.l 


mask5 * 


14.4 


CDFS.l 


mask6* 


12.0 


CDFS.l 


mask7 


19.2 


CDFS.2 


mask8 


19.2 


CDFS.2 


mask9 


16.8 


CDFS.2 


mask 10 


9.6 


CDFS.l 


maskl 1 


13.2 


CDFS.l 


mask 12 


15.6 


CDFS.l 


mask 13 


9.6 


FIRES 


mask 14 


16.8 


FIRES 


mask 15 


13.2 


FIRES 


mask 16 


10.8 


FIRES 


mask21 


15.6 


CDFS.2 


mask22 


18.0 


CDFS.2 


mask23 


15.6 


CDFS.2 



geometrical constraints, once all the highest priority targets were 
assigned, the mask was filled with lower priority targets (typically 
more uncertain photometric redshifts or K magnitudes outside our 
main sample). The total weight for all targets allocated in the mask 
was then calculated and a Monte-Carlo technique used to design 
many realisations, keeping the version with the highest weight. 
Multiple masks at the same pointing position were then designed, 
giving highest priority to those galaxies not already observed in a 
previous mask (but repeats of some slits occur and are useful for 
checking purposes, discussed later). The exact details of the object 
prioritising are unimportant as a) the completeness is calculated a 
posteriori in £|3 - 3 1 and b) eventually enough targets are observed to 
ensure high completeness for all galaxies within this K magnitude 
range after ~6 masks. The main reason for this prioritisation is to 
ensure a relatively high completeness after only a few (~3) masks 
have been observed. 

We observed each mask for typically 20 N&S cycles of 60s 
exposures at each of the A and B positions, resulting in an expo- 
sure time of 40 minutes for each image. We typically took 6 expo- 
sures resulting in a total integration time of four hours. A log of the 
observations is given in Table Q] The seeing ranged from 0.6" to 
1.0" and was typically around 0.8". It should be noted that the few 
frames with significantly shorter exposure times are lower prior- 
ity masks and so the overall loss in completeness if these masks are 
slightly shallower is minimal. The final numbers of slits we are able 
to place on ROLES' targets (i.e., matching our A"-band selection 
criterion and photo-z weighting) total 1849 for CDFS and 533 for 
FIRES. These totals refer to the number of spectra we might expect 
to obtain for unique objects (i.e., repeated slits are only counted 
once) after discounting 'filler targets' from other programmes and 
the few objects (<5 per mask) which resulted in corrupt spectra 
due to poorly-cut slits, slit collision or contamination by scattered 
light, etc. 



2.3 Spectroscopic data reduction 

2.3.1 Pre-processing 

In addition to standard spectral CCD data pre-processing (overscan 
subtraction, bad pixel masking etc.) a number of non-standard steps 
are necessary for our data, which we detail here. We used several 
routines from the Carnegie Observatories System for Multi-Object 
Spectroscopy 2 (COSMOS20), but the majority were applied using 
custom written routines in IDL. Examples of three representative 
slits from one mask are shown at various stages of the reduction in 
Fig.[3] In all cases the same portion of the frame is shown, which is 
about 100 pixels (~270A) long in the spectral direction. Panel "a" 
shows a single, raw exposure prior to sky subtraction. 

Firstly the COSMOS2 STITCH routine is used to combine the 
images from the two separate LDSS3 amplifiers into a single frame 
and convert the counts in the image into electrons. Overscan sub- 
traction and bad pixel masking were then performed (the bad pixels 
being identified by significantly deviant pixels seen in the ratio of 
two dispersed sky flats). 



2.3.2 Charge trap and bad pixel masking 

The repeated charge shuffling involved with N&S observations can 
cause artefacts which appear as short streaks along certain columns 
of the CCD. This is due to localised defects ('charge traps') in the 
CCD which cause the charge to become smeared as it passes over 
them. Charge traps are identified by taking dark frames using the 
same N&S strategy as for the science data. This leads to a frame 
containing columns of bright pixels along the shuffle direction due 
to charge trap defects. These are identified by median smoothing 
the image using a 1 x 3 pixel boxcar (so that the long axis is per- 
pendicular to the shuffle direction) and taking the ratio of the origi- 
nal N&S dark to the median-smoothed version. Charge traps show 
up as significantly deviant pixels and are flagged and added to the 
original bad pixel mask. After using this charge trap mask only a 
very small number of obvious charge traps still remain and we dis- 
cuss how these are dealt with in ^2.41 The same charge trap mask 
is used for data from all observing runs, since the pattern appeared 
stable and the few residual charge traps are removed by hand, as 
described below. 



2.3.3 Registration and wavelength calibration 

The next step is to calculate the transformation between each in- 
dividual exposure, to register all exposures of the same mask to 
a common frame. There are two consequences leading to offsets 
caused by flexure within the telescope and/or instrument (or not 
perfectly reproducible positioning of the slit mask within the in- 
strument). The first means that the images of the slits may not 
fall on the same pixels each time. This is especially the case when 
data were taken on different observing runs. To deal with this, we 
use the mapping routines in COSMOS (MAP-SPECTRA, ADJUST- 
MAP, etc.) to accurately identify the locations of a set of reference 
sky lines in every single exposure. By taking the first frame as a ref- 
erence and requiring the positions of the sky lines to match in each 
subsequent exposure, we can calculate the linear shift and rotation 



2 see http://users.ociw.edu/oemler/COSMOS2/COSMOS2.html 
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Figure 3. Examples of the various steps in the data reduction, showing a region containing three representative slits. The six panels each show a region of 
100 pixels (270A) in the spectral (horizontal) direction by 100 pixels in the spatial (vertical) direction. Panel a): the initial individual frame. Panel b): the 
same frame after sky-subtraction (by shifting a copy of the image by 16 pixels and subtracting) Panel c): after median combination of several exposures with 
cosmic ray rejection. Panel d): coaddition of the positive and negative (A'-B') spectra after subtracting the two different nod positions Panel e): the noise 
frame calculated as described in the text. The highest levels of noise correspond to regions underneath bright sky lines. Panel f): The final signal-to-noise map. 
Black corresponds to a level of 3. Oct and white to 7.0<r. Two peaks corresponding to emission lines clearly stand out. Tracing them back to earlier stages of the 
reduction, the line in the lowest slit can be seen in an individual exposure (panel b), but the line in the central slit only becomes vsisible after stacking several 
frames (panel c). The dipole signature of a genuine emission line is clearly visible in the two individual nod position spectra (A' and B'). 



(and occasionally a small scale change) to register the imagefl. The 
second possible offset involves the movement of the galaxy within 
the slit and will be dealt with later. Up to this point, we have calcu- 
lated how to align the images of the sky spectra. 

ADJUST-MAP in COSMOS does a good job of accurately cen- 
troiding on the sky lines in the dispersion direction, however signif- 
icant offsets are seen from the centre of the slit in the spatial direc- 
tion. This may be because the routine is not optimised to work with 
N&S data. Since accurate slit location is essential (as described 
later) we use our own slit tracing routine to overcome this. Briefly, 
this code takes many ID profiles through each spectrum in the spa- 
tial direction and cross-correlates this against a model of a N&S 
spectrum (approximately two top-hats of width 16 pixels, the slit- 
width/shuffle distance, side-by-side) to find the optimum offset. In 
the tracing and rectification, slit curvature and tilt are neglected. 
Since the spectra are short (370 pixels), the change in y coordi- 
nate over this length is typically one pixel or less, so spectra are 
extracted in rectangular boxes aligned to the CCD axes, for sim- 
plicity. 

Before applying this transformation and resampling the spec- 
tra, it is necessary to perform the sky subtraction as described be- 
low. If this is not done, a small rotation can result in strong sky- 



3 Since LDSS3 is oversampled, we use nearest-neighbour resampling for 
all transformations, to avoid correlating the noise between neighbouring 
pixels and hence compromising the noise estimate in each pixel. 



subtraction residuals due to the sharp edges of bright sky lines be- 
ing under/over-subtracted as they pass into a neighbouring pixel. 

As a by-product of calculating this transformation, maps are 
constructed to convert global x and y CCD coordinates to wave- 
length and slit spatial position (specifically slit centre). These maps 
will later be used to identify the wavelength of a detected emission 
line, transforming from its CCD coordinates. 

2.3.4 Sky-subtraction 

For each individual exposure (see Fig. [3] panel a), N&S sky sub- 
traction is performed by shifting the image by the shuffle distance 
(16 pixels) in the spatial direction and subtracting it from itself 
(Fig. [3] panel b). If we consider A-B, in the notation introduced 
in the previous section, every slit will now contain a positive image 
of galaxy spectrum A' (where the prime is used to denote a sky- 
subtracted spectrum) and a negative image of spectrum B' of the 
same galaxy, both with the sky removed. 

Once the sky subtraction has been performed the spectra are 
resampled to a common frame, as described above. Using the 
boundaries identified by the slit-tracing program, all pixels outside 
the spectra are flagged as bad and their values omitted from further 
analysis. This is crucial to avoid contamination of the edges of re- 
gions of good data with the bright sky spectra at the edge of the 
slit. This also ensures that the correct exposure time is propagated 
(by tracking the number of good exposures in each pixel after the 
next alignment step). The next step is to align the positions of the 
galaxies within the slits. Although the positions of the slits have 
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been registered, the position of the galaxy can shift by up to a cou- 
ple of pixels due to alignment drift of the instrument. To register 
the positions of the object spectra, several bright nebular emission 
lines in the galaxy spectra are identified in each mask. The cen- 
troids of these lines are measured in each individual exposure and 
offsets calculated^ For all masks, the shifts are found to only be 
significant in the y (spatial) direction (confirming the accuracy of 
the wavelength calibration) and the spectra are then resampled to 
the nearest integer pixel shift. 

Next, all the individual exposures are median-combined with 
5a outlier-rejection (where a is the noise estimate, described 
below) to remove cosmic rays (Fig. [3] panel c). (Experimenta- 
tion showed that median combination gave a slightly cleaner co- 
addition than using the mean.) In order to be left with a positive, 
summed image of the galaxies, the sky-subtracted spectra can sim- 
ply be shifted by their separation, 6 pixels (which is equal to the 
nod distance, 1.2"), and subtracted, A'-B' (Fig. [3] panel d). Notice 
that in Fig. [3] panel d, pixels outside the region of the combined 
science spectra are set to bad values (appearing blank in the plot). 
Prior to object-finding, the extracted 2D spectra are shrunk by one 
pixel in the spatial direction (i.e., above and below) on either side of 
the slit. The nearest neighbour resampling technique means that all 
transformations are accurate at the one pixel level, and so this pre- 
vents rounding error. However, it can be seen that in the examples 
shown, bright sky lines are bleeding into the lower portion of the 
slit (appearing black in the image, since they have been subtracted 
in the previous step). Shrinking the extraction window helps to re- 
duce the impact of such artefacts. However, since the spatial extent 
of the emission line is very close to the width of this window, care 
must be taken to compensate for flux which may fall outside this 
window (the top edge in this case), discussed in £12.5 . 1 1 

This then leaves us with the final, reduced 2D spectr^fl from 
which we will make measurements such as line flux. However, in 
order to optimally detect emission lines in the first place, we will 
further process the images to produce optimally filtered versions. 
We also require an estimate of the noise in each pixel of our detec- 
tion images to assess the significance of our detections, which are 
described below. 



4 Since only a small number of objects per mask (<6) possessed emission 
lines bright enough to be found visually in an individual exposure, we were 
unable to fit for any transformation other than a small shift. It is possible 
that the galaxy positions may be rotated between exposures, for example. To 
assess the accuracy of this transformation, we calculated shifts from three 
bright lines and looked at residual shifts in other bright lines in the same 
exposure. Typically these were centred to better than 2 pixels. 

5 One step which has been intentionally omitted is flatfielding. Tests run- 
ning the data through our whole analysis show that attempting to flatfield 
makes negligible difference (the N&S methodology will also naturally aver- 
age over some pixel-to-pixel sensitivity variations) and possibly marginally 
increases the noise. The main limitation to the accuracy of the flatfielding is 
the accuracy with which the flatfields can be registered to the science data. 
The flatfields cannot be registered in the same manner as the science data 
since the latter rely on the 2D position of night sky lines. The flatfields may 
only be wavelength calibrated from arc lamps, taken periodically through 
the night, and by the positions of the slit edges to locate the spatial off- 
set. This can lead to small systematic offsets between the flatfield and the 
science data which may lead to increased noise when the (offset) flat is 
applied. 



2.3.5 Noise estimation 

By following a similar procedure to the sky-subtraction just de- 
scribed (i.e. shifting a copy of the image by the shuffle distance), 
but adding instead of subtracting the frames, an image is produced 
which contains a spectrum of the sky instead of the galaxy at the 
positions of the galaxies (A and B) in each slit. By again median- 
combining the individual exposures a frame is created with the best 
estimate of the value of the sky at each pixel. This is used to es- 
timate the noise, N inaiv , in each pixel of each nod position in an 
individual exposure. The noise in the i, j th pixel is given by 



= V\ < sky > \ij + 2R 2 , 



(1) 



where < sky > is the median-combined sky frame, R is the 
read noise of the detector (3.4 electrons), and the factor of 2 ac- 
counts for the fact that two readouts have been combined in the 
sky-subtraction/addition by using the two shuffle positions. 

After stacking n frames exposures, the noise in the i, j pixel, 



N C om,ij ^indiv.ij j yj Ttfrar. 



(2) 



Once the two nod positions in the combined frame produce an 
image containing the sum of the objects, the noise in the i, j pixel 
of the final combined science image is then 



IV. . — I ]\T2 I AJ2 

l J \j com j^r ,ij com^t^iji 



(3) 



which is just the quadrature sum of the A' and B' position noise 
values (Fig. [5] panel e). 



2.4 Line Detection 

2.4.1 Optimal filtering 

We wish to search the 2D spectra objectively for significant fea- 
tures cc^tsistoit_wdA emission lines from the target galaxies (see 
also iGlazebrook et alj|2004bl for a similar application). Since the 
optimal detection kernel is the profile of the line itself, we start 
by finding several obvious bright lines, easily visible by eye, and 
measuring their profiles. These are typically well-approximated by 
elliptical Gaussians with FWHM of 5 pixels in the spectral direc- 
tion and 3.5 pixels in the spatial direction. The kernel is normalised 
by setting its total to unity to conserve flux. 

The science images are convolved with this kernel, producing 
a signal frame, S (Fig. [3] panel f). This convolution changes the 



properties of the noise in each pixel from Nij to N c , 



where 



N conVyij = JN? j g> k 2 , 



(4) 



where k is the emission line-shaped kernel described above. 
N conVt ij is thus the quadrature sum of the contributions of each 
pixel in the kernel to the noise in the i,j pixel. 

Now, for the i, j th pixel, the significance, a, is given by 



5 •.! 



(5) 



J N 

where by is the local background (continuum) estimated from the 
mean of all pixels in two ID side-bands, each 40 pixels wide, rang- 
ing from i + 10 to i + 50 pixels and i — 10 to i — 50 (or whatever 
number of pixels from these regions falls on pixels of the science 
spectrum flagged as good). It is appropriate to only measure the 
continuum in a ID box (i.e. a line at the position of the j th pixel) 
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in this way, since the image has already been convolved with a ker- 
nel in the spatial (j) direction and thus contains an estimate using 
values from multiple pixels in this direction. This simple approach 
to continuum estimation has the drawback that it can be severely bi- 
ased high when a bright emission line enters one of the sidebands. 
To alleviate this problem, we perform one iteration, re-estimating 
the continuum in the same way after masking out allpixels above 
a high significance threshold (6a) from the first pass[j 

The above is an exact expression considering sky noise and 
read noise. However, in some cases we cannot neglect the contri- 
bution to the noise from the Poisson noise in the continuum (even 
for these faint objects). We adopted the following approach as an 
operational way to consider the effects of continuum noise. 

The maximum contribution from the continuum is taken as 
the maximum within one pixel above or below the position of the 
candidate emission line: 



bmax — max{bij~i,bij,bij + i} 



(6) 



(to allow for the effect of rounding error in locating the position of 
the continuum). Then, the noise in the kernel-smoothed image is 
estimated as 



bnoise — \/ bij ® k 2 3 



(7) 



where the bij term here is from the Poisson noise squared. The final 
significance estimate is then modified to 



Si 



V conv 



. +b 2 . 

tj 1 noise 



(8) 



for sources with significant flux in the continuum. Significant flux 
can be as little as ~5-10 electrons, since other sources of noise, 
even bright sky lines, have been beaten down by repeated observa- 
tions. 

We calculate the significance of each pixel in the science spec- 
tra and locate the most significant peaks. Such a significance map 
is shown in Fig. [3] panel f. Significant pixels are grouped together 
using a standard n-connected pixels algorithm (joining all neigh- 
bouring pixels above the minimum significance threshold into a 
single group), identifying the peak by the single most significant 
pixel. In addition, detections are rejected where the peak pixel lies 
at the edge of the spectral extraction box, as these are always due 
to edge effects from bright sky lines. 

2.4.2 Catalogue cleaning 

The majority of artefacts present in the emission line catalogue 
come from one of two sources. There are: bright sky lines leak- 
ing into the edge of the slit extraction box, and broad 'lines' caused 
by weak continuum emission which is not well sampled by our 
simple sliding box continuum filter. These all needed to be flagged 
interactively. 

Fig.[4]shows examples of convolved and un-convolved 2D and 
ID spectra for a range of line significances. Such plots are visually 
inspected, sorted by significance, for all objects. Only detections 
which are due to obvious artefacts (edge effects from sky lines, 
missed charge traps, scattered light from bright neighbouring spec- 
tra, underestimated continuum) are rejected. We make no additional 

e This threshold was determined by experimentation. It was found that 
emission lines with an initial significance <6ct had negligible impact on 
the estimate of their surrounding continuum, and thus this threshold is a 
conservative threshold for re-estimating the continuum level. 



cuts on a subjective basis, such as line shape, for example. A rea- 
sonably clear cut-off in the reliability of lines is seen at ~ 4a in 
our units using this technique, but this will be re-assessed quantita- 
tively in the next section. Above this limit, the shapes of the lines 
in all the plots appeared broadly consistent with each other. After 
cleaning the catalogue, we are left with emission line detections for 
246 unique objects in CDFS and 1 19 in FIRES, from the 1849 and 
533 unique objects (CDFS and FIRES, respectively) on which we 
initially placed slits. There are several tests which can be performed 
to assess the reliability of our emission line catalogue. 



2.4.3 Tests of emission line catalogue 

The most basic quality check is to compare the properties of emis- 
sion lines detected in the catalogue with repeat observations of the 
same galaxy observed in a different mask. In this comparison we 
select galaxies from our emission line catalogue and identify which 
of these have been observed in more than one mask. 

There are 44 such repeat slits between these masks (observed 
a total of 100 times - some objects are observed on several differ- 
ent masks). Each detection is tested in turn to see if it has been 
reproduced in a repeat observation of another mask. Each detection 
is considered individually (some z~0.5 galaxies have as many as 
three genuine emission lines in each slit). If it is independently de- 
tected on every mask, we consider this successfully reproduced. If 
it is not found on every mask, it is added to a list of non-reproduced 
objects. We note that this is conservative, as an object indepen- 
dently discovered in two out of three masks would be classed as 
spurious in this scheme. Also, there may be some overcounting in 
this method as it is possible that if a galaxy is observed in two 
masks and both produced discordant detections, both would be 
counted as spurious when in fact the more significant line may be 
genuine for the following reason. Due to the slightly different flux 
limits of the different masks, it is possible that a line observed in 
one mask may be just below the detection threshold in a shallower 
mask. We make no attempt here to correct for the different depths, 
and thus a line overlapping with a slightly shallower mask may be 
undetected in the shallower mask and thus incorrectly flagged as 
spurious. Regardless, using these criteria, the objects successfully 
recovered are plotted in Fig.[5](solid line) and those not recovered 
(dashed line). The cumulative fraction (main panel) shows that 95% 
of all false detection^ occur below 4.0cr. We reiterate that this is 
likely an upper limit, for the reasons just discussed. 

This internal comparison also allows an estimate of the accu- 
racy of the wavelength calibration. The difference in wavelength 
between lines successfully recovered is well described by a Gaus- 
sian of width ct=1.4A, which is ~0.5 CCD pixels. This corresponds 
to an observed frame velocity difference of ~100 km s _1 for 
ROLES galaxies or ~50 km s _1 rest frame at z~l. 

A second internal test of the catalogue is possible for those 
galaxies with multiple emission lines. We can attempt to measure 
redshifts from multiple lines and check that the wavelengths of the 
lines detected are consistent with a single redshift for the object. 
This is done by considering all possible combinations of common 
galaxy emission lines, compared with the observed wavelengths of 
our detected lines. Plots of these test redshifts for each galaxy are 
visually inspected and the most likely combination assessed. In al- 
most all cases there is either one or zero match. For a small number 
of cases, implausible combinations of lines are rejected (such as 



The line detection algorithm stores all detections down to 3. Oct. 
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Figure 4. Examples of 1- and 2-D spectra, ranked in order of decreasing significance. Each panel shows the smoothed ID spectrum (thick black line), centred 
on the detection, as the main plot. The solid line bounded by two dashed lines shows the estimate of the continuum and the lcr noise limits around this, 
respectively, and the dotted line indicates the nominal 4cr threshold. The plot is annotated with the significance of the detection. The two insets below show 
the smoothed (upper inset) and unsmoothed (lower inset) 2D spectrum. The wavelength scale is shown below. 




Figure 5. Reproducibility of lines from repeated observations of the same 
galaxies in mutliple masks. Solid lines show those lines successfully re- 
produced in multiple masks, dashed lines show objects which weren't suc- 
cessfully independently found on every repeated mask. See text for further 
details. Main panel shows the cumulative fraction of such detections, indi- 
cating that 95% of 'spurious' identifications occur at significance <4.0o\ 
Inset shows a histogram of the number of detections, zoomed in to the low 
significance end of the plot. 



[N II], [S II], but no Ha). Other than this, no astrophysical consid- 
erations of likely line ratios are used to make the decision. Where 
no plausible combination is found, the most significant line is re- 
tained as a single line detection and all lower significance lines in 
that slit are rejected. This may not be a completely representative 
number for the whole survey as some slits contribute several lines 
to the spurious count by this measure. To reduce the impact of the 
few slits with many false positives, we only count a maximum of 
one false positive per slit. This is likely a better estimate of the true 
rate in our survey, since most slits have at most one detection. In 
a similar manner to Fig. [5] Fig. [6] shows the likely spurious detec- 
tion rate arising from non-physical combinations of lines. Note that 
90% of spurious lines by this measure occur below a significance 
of 4.0cr. 

75 objects result in plausible multiple-line redshifts. Of these, 
61 belong to z<< 1 galaxies (mostly z~0.5); 14 show [O II] at z~l 
plus at least one other emission line. 

As a further test for these multiple line objects (secure red- 
shifts) we can compare the redshifts from our data with others avail- 
able from the literature. Where possible, public spectr oscopic red- 
shifts are used. For CDFS we use the compilation from Wu vts et"al] 
d2008h and only consider secure redshifts (their quality flag 1.0). 
This CDFS sample comes from a variety of surveys with very dif- 
ferent selection criteria and thus caution must be used when con- 
sidering how representative of the whole of ROLES are any statis- 
tics derived from a comparison with these redshifts. For the FIRES 



ROLES II: method and SFRD 9 




Figure 6. Using the combination of multiple lines within each slit to as- 
sess the genuine/spurious nature of detections. Solid lines show slits where 
the combination of lines yield a self-consistent redshift, dashed lines show 
those which do not. Main panel shows the cumulative fraction of such de- 
tections, indicating that 90% of 'spurious' identifications identified in this 
way occur at significance <4.0<r. Inset shows a histogram of the number of 
detections, zoomed in to the low significance end of the plot. See text for 
details. 



field, we use spectroscopy from lTran et alj |2007) (kindly provided 
by K.-V. Tran) primarily designed to study the z=0.83 cluster, 
MS 1054-03, below the redshift limit of ROLES. Again, due to the 
very different selection criteria from ROLES' galaxies, care should 
be exercised in drawing comparisons. A plot of the ROLES red- 
shifts against these others is shown in Fig. [7] Points with error bars 
denote photometric redshifts and red squares denote secure pub- 
lic spectroscopic redshifts. Three secure public spectroscopic red- 
shifts disagree with our multiple line spectroscopic redshifts. For 
two of these, we assign a redshift of z~0.5 versus z~l for the pub- 
lic redshifts. For these two cases, a second spurious line is seen in 
ROLES data which makes a genuine [O II] emission line appear to 
be H/3 or [O III]. For the third discrepant point, the public redshift 
puts [O II] just outside our redshift window, meaning that we have 
assigned two spurious lines to [O II] and [Ne III] when we should 
actually have recorded no emission lines. These multiple spurious 
lines coinciding with the correct spacings for genuine lines are rare 
in ROLES as a whole0 For the remaining secure public spectro- 
scopic redshifts, we agree with the public redshift. Similarly, the 
vast majority of photometric redshifts agree well with our multi- 
ple line spectroscopic redshifts. A small number of outliers exist, 
which may be due either to serendipitous combinations of spurious 
lines, or catastrophic failures of the photometric redshift. We note 
that plotting a single best fit photometric redshift and confidence 
interval is not quite the same as the method we use to interpret the 
redshift of an emission line (as described below) and, in practice, 
our weighting scheme will be more robust than this simple compar- 
ison. 



8 For the few rare instances where a secure public spectroscopic redshift 
exists which disagrees with a ROLES' redshift, we favour the former value, 
since the alternate spectroscopy likely has multiple emission lines (due to 
the longer wavelength coverage than our band-limited data) or other fea- 
tures which make the other redshift more secure than ours. 
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Figure 7. Comparison of secure spectroscopic redshifts from ROLES 
galaxies possessing multiple emission lines with independent redshift mea- 
surements from the literature. Filled circles with error bars denote photo- 
metric redshifts from the FIREWORKS and FIRES catalogues; red squares 
denote secure spectroscopic redshifts from the sources discussed in the text; 
blue squares show lower confidence (^ 0.5) public spectroscopic redshifts. 
The vast majority of public redshifts agree well with our multiple line red- 
shifts, and we discuss in the text the reasons for the small number of outliers 
(which are likely more common in this comparison subsample using mul- 
tiple lines than for ROLES as a whole, where the typical object exhibits a 
single line). 



2.5 Line identification 

For the majority of these emission line objects, only a single line 
is detected and thus the spectrum on its own gives an ambiguous 
redshift. Thus we must make use of additional data to determine a 
redshift. For a subsample of these objects, spectroscopic redshifts 
are available from public spectr osco py. We aga i n use the compi- 
lations from IWuvts etalJ d2008h and iTran et a U ( |2007|) . consider- 
ing only secure redshifts. For most objects, public spectroscopic 
redshifts are not available (only 24 ROLES' galaxies already have 
spectroscopic redshifts from public spectroscopy, which is not sur- 
prising given their faintness) and we must use photometric red- 
shifts in order to determine the most likely identity of our emis- 
sion line. Since our survey began, more accurate public photomet- 
ric redshifts for th e CDFS have become available (FIREWORKS, 
IWuvts et al.l2008h . They utilise many more filters and cover a wider 
wavelength range than those used in the Mobasher & Dahlen cata- 
logues and hereafter we adopt these. For consistency with the orig- 
inal /('-band selection, we retain the original i^-band total mag- 
nitudes from the Mobasher & Dahlen cata logue^ For the FIRES 
field, the photometric redshifts derived in Fors ter Schreiber et al.l 
( 2006) have been used throughout. The most probable identities of 



9 The agreement between total magnitudes is generally very good since 
they are based on the same original images, but with a slightly different 
reduction and photometry procedure. However, a small number of sources 
from the Mobasher & Dahlen catalogue are not present in the Wuyts et al. 
catalogue. This is presumably due to the former being fi-band selected and 
the latter being A'-band selected, and also the A"-band surface brightness 
limit imposed by the latter and/or slightly different masking/deblending pa- 
rameters used in the object detection. 
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Figure 8. Photometric redshift probability distribution functions for three 
example galaxies. The probability is in arbitrary units. Shaded regions 
show the redshift windows for [O II] (rightmost region, darker shading); 
H/3/[0 III] (centre shaded region); and Ha (lowest redshift region). The 
top panel shows a galaxy where the detected emission line would be con- 
sidered [Oil], the centre panel one which would be considered definitely 
not [O II], and the lower panel shows a galaxy which falls between the two 
cases and has a probability Pon ~ 0.3. 

strong emission lines in our survey are [O II]A3727 (0.889 < z < 
1.149), HaA6563 (0.073 < z < 0.220) or one of H/3A4861 or 
[OIII]AA4959,5007 (0.406 < z < 0.648). Since we only care 
whether or not the line is [O II] we consider the likelihood that 
it is [O II] versus one of the other lines. In order to do this, the 
probability distribution functions (PDFs) of the photometric red- 
shifts for both fields are used, kindly provided to us by S. Wuyts. 
The photo-z PDF is integrated over the redshift range for which 
[O II] could be detected and this is compared with the PDF inte- 
grated over the three different redshift windows just described. The 
ratio then gives the probability that the line is [O II] (Pon). If the 
probability is Pon <0.1 it is set to zero (i.e. assume it is not [O II]) 
and for Pon >0.9 it is set to unity (assume it is definitely [O II]). 
The majority of our detections fall into one of these categories. 

Of the 246 detected emission line objects in CDFS and 1 19 in 
FIRES, only 17 have 0.K Pon <0.9 in CDFS and 40 in FIRES. 
195 have Pon >0.9 in CDFS and 50 in FIRES. The higher frac- 
tion of ambiguous objects in FIRES relative to CDFS is due to the 
slightly less well-constrained photo-z PDFs of the former due to 
the different photometric filters available for each field (see §3.1). 
For the few cases where 0.1 < Pon < 0.9 this weighting is prop- 
agated through the analysis such that these galaxies contribute a 
fraction Pon of their properties to the measurement under consid- 
eration (i.e. a galaxy with Pon = 0.5 would contribute half its 
[O II] flux to a measurement of the total [O II] flux). Fig[8]shows 
examples of PDFs for the three different cases. 

2.5.1 Flux measurement 

Line fluxes from the 2D science frames are measured in a 5 x 5 
pixel box around the centre of each emission line. For a few detec- 
tions, one row of pixels in this box may fall outside our masked re- 
gion (i.e., onto bad pixels). In these cases the profile is corrected by 
mirroring pixels on the opposite side of the 5x5 box into the bad 
region. Flux errors are determined by applying the same method 
to the noise frame (without convolution). It is important to distin- 
guish between the detection errors/significances, computed from 
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Figure 9. Comparison of line fluxes between IVanzella et alj j2008l) 
(FORS2) and ROLES for all objects in common. Numbers indicate ROLES 
mask number. The two biggest outliers lie in the O2 A-band region - the 
ROLES fluxes are corrected for telluric absorption, whereas the FORS2 
fluxes are not. 



the optimally-filtered noise frame, and the flux errors, computed 
from the original (i.e. prior to convolution) noise frame resulting 
from eqn.[3] 

To calibrate these fluxes, the ESO spectrophotometric stan- 
dard star LTT3864 was observed with the same instrumental setup, 
through one of the slits of the science masks. To compensate for 
absorption by telluric features, we isolate the O2 A-band feature 
around 7600A in our flux standard and measure the ratio of the 
observed flux to that expected from the smooth fit. 

The zeropoint of our flux scale is verified using measurements 
of science targets overlapping with the extensive spectroscopy in 
the CDFS. We use ESO public ID spectnQfl observed with FORS2 
in the CDFS JVanzella et al. 2008). Objects in common with our 
survey are identified and line fluxes are measured using a similar 
method to that used in ROLES (but on the ID spectra). The flux 
calibration of the ESO public data has been extensively checked 
against multi-band ACS photometry for a large sample of objects 
and there does not appear to be anything unusual about the ob- 
jects in our ove rlapping subsample (E. Vanzella, priv. comm.). The 
IVanzella et al. I d2008h spectroscopy is calibrated to reproduce the 
total flux in a point source. We have examined the sizes of our 
ROLES galaxies in ACS i-band images and find that most will have 
sizes very close to a point source when o bserved with our typi- 
cal ground-based seeing (~0.8") and thus IVanzella et alj ( 2008) 
should give us approximately total fluxes for our objects. We plot 
the flux in our spectroscopy against the public FORS2 spectroscopy 
in Fig. [9] We find excellent agreement between the two samples, 
with the scatter dominated by measurement uncertainties. Error 
bars have been omitted for clarity, but the size of our errors can 
be gauged from Fig.[l0]and the FORS2 errors are comparable. The 
numbers indicate the ROLES mask on which our measurement was 
made. There are no obvious systematic offsets with mask num- 
ber, indicating that this single zeropoint calibration should apply 
equally well to our whole survey. 

The repeatability of our flux measurements is shown in Fig.llOl 
from repeat observations as described in ^2.4.31 The agreement be- 
tween independent measurements of the flux in different masks is 



see: http://www.eso.org/science/goods/ 
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Figure 10. Flux difference between repeated measurements of emission 
lines on different masks. Inset panel shows the distribution of flux difference 
divided by the combined flux errors. The overplotted curve is a Gaussian of 
width 0.9<7 indicating that our estimates of the flux errors are reasonable. 
See text for discussion. 



excellent. In the inset panel, we compare the errors on each mea- 
surement with the combined errors (sum in quadrature) of these 
repeated observations. We find that the resulting distribution is ap- 
proximately Gaussian, with a width of 0.9a indicating that our es- 
timates of the noise are reasonable. 



3 ANALYSIS 
3.1 Stellar masses 

Stellar massed] are obtained by fitting the extensive 
multiwavelength photometry at the spectroscopic red- 
shift to a grid of stellar popul at ion m odels (using PE- 
GASE.2, Fioc & Rocca-Volmerangel 1 19971) as described in 
iGlazebrook et all b004al) . For the CDFS, this consists of 
C^38, B435, B, V, Veoe, R, i775,I, Z850, J, H K a , [3.6fi m], 
[4.5/um], and [8.0/xm] photometrjE] (see IWuvts etafl 120081 for 
details); and for FIRES t his consists of U, B, V, Vpoe, Isn, Js, H, 
and K s photometry (see iForster S chreiber et al. 2006 for details). 
The FIRES c atalogue was co r rected for Galactic extinction using 
the maps of ISchlegel et all dl998h . For CDFS the extinction, 
even in the blue, is negligible. The matched apertures used for 
measuring colours (see the cited papers for details) are used 
in the SED-fitting. Fitted quantities related to total luminosity 
in the stellar population models, such as stellar mass and SFR, 
must therefore be scaled from the aperture measurements to 
total light measurements. This is achieved by multiplying by 
lQ-0-i(K tot -K ap ) wnere Kap and KtQt ^ the ^'_band aperture 

magnitude and the total magnitude, respectively. 

We note that the colour distributions for the CDFS and FIRES 
fields are significantly different (both for our spectroscopic sample 
and for the parent catalogues, although this is consistent with the 
expectations of cosmic variance for our narrow redshift slice) in the 



11 Stellar masses throughout this paper refer to the present mass in stars, 
not including remnants 

12 In paper I we used a subsample of these filters, not including the Spitzer 
data. As we noted in that paper, omitting the Spitzer data makes negligible 
difference to the fitted masses at these redshifts and thus the absence of MIR 
data for FIRES does not compromise the stellar masses for this field. 



sense that CDFS galaxies are bluer. This means that for a given K- 
band magnitude, CDFS probes to lower stellar masses than FIRES 
by around 0.2dex. 



3.2 Nominal L([OII]) Conversion 

We use the same conversion of L([ OII]) to SFR as J05 which in 
turn was based on the calibration of lKennicuttl il998l) . We convert 
L([OII]) to L(Ha) assuming that {[Ott\IYLa) obs = 0.5. We then cor- 
rect for extinction in Ha assuming average extinction of Aij Q =l. 
Thus the conversion is: 



SFR{M Q yr- 



10° 



0.5 



7.9 x 10" 
1.82 



-Lipimergs" 1 ) (9) 



where the factor of 1 .82 accounts for the conversion from a Salpeter 
IMF JSalpeteilll955b to that of BG03. 



3.3 Spectroscopic completeness 

In more traditional spectroscopic surveys, it may be more natu- 
ral to split the completeness into 'targeting completeness' (i.e. the 
fraction of galaxies of interest on which slits are placed) and 'red- 
shift success rate' (the fraction of objects with slits yielding red- 
shift). Since we use photo-z selection and implicitly assume we ob- 
tain redshifts for 100% of objects within our redshift window with 
[O II] emission above our flux limit, we simply combine these into 
a single 'completeness' term. To calculate the completeness of our 
spectroscopy, we again use the FIREWORKS and FIRES photo- 
metric redshift PDFs and public spectroscopy. For every galaxy in 
our A"-selected sample, in the area of our spectroscopic pointing, 
we sum their PDFs. This summation gives us the total redshift dis- 
tribution of all galaxies in the sample (upper curve of Fig. II It. The 
number of galaxies in the k th magnitude bin in the redshift range, 
zo to zi is then 



N k = 



Pk(z)dz 



(10) 



where P(z) is the redshift probability density function given 
by the photo-z PDF.0 

We then repeat this process for just the galaxies on which we 
placed slits. This gives the lower, thicker, green curve in Fig. QT] 
Since we are only interested in a specific redshift range (the shaded 
region in the figure), the completeness is given by the ratio of the 
integrals of the PDFs over this redshift range, i.e. the ratio of the 
area under the upper curve to the lower curve in this redshift range, 
i.e. 



Wk 



(11) 



phot 



where Wk is the completeness in the k th magnitude bin, N s m b 
is the number of galaxies on which we placed slits and N v hot is the 
number in the entire photometric catalogue, in the same field, and 
N are calculated as given in eqn.1101 

Our overall survey completeness is ~0.85 for CDFS and ~0.7 
for FIRES. In practice, the completeness is a function of A'-band 
magnitude, so we repeat this procedure in 0.25 magnitude bins. 
For CDFS, our selection has removed this dependence, making it 

13 Replacing the photo-z P(z) with a delta function at the spectroscopic 
redshift, where redshifts are available from public spectroscopy, makes neg- 
ligible difference to the resulting completeness. 
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Figure 11. Redshift distributions used to calculate the completeness in 
FIRES (upper panel) CFDS (lower panel). These are constructed by sim- 
ply summing the FIREWORKS photometric redshift P(z)'s for individual 
galaxies (or delta functions for galaxies with spectroscopic redshifts). Plot 
shows distributions for all 22.5 < K ^ 24.0 galaxies in our survey area 
(upper black line) and those on which we placed slits (lower green line). 
Vertical lines indicate the redshift limits of our survey. The completeness is 
given by the ratio of the integral of the upper line to the lower line within 
this window. See text for details. 



FIRES 



0.41 

22.5 



Figure 12. Spectroscopic completeness, calculated as described in the text, 
as a function of i^-band magnitude for the two ROLES fields. Note that the 
completeness axis shows down to 40%, not 0%. 



nearly independent of magnitude, but for FIRES there is still a weak 
magnitude dependence, the completeness dropping to 0.50 in the 
faintest binl 14 lThe results of this are plotted in Fig.l 121 



3.4 Survey volume 

In order to calculate the SFRD (i j4.2l ( we require an estimate of the 
maximum volume, Vmax, from which each galaxy observed in our 
survey is drawn. For each galaxy, the value of Vmax is determined 
by two limits: an [O II] line flux limit and a 7\T-band flux limit. 



14 The lower overall completeness for FIRES relative to CDFS for the 
same observing strategy is due to the presence of the cluster MS 1054-03 
at z=0.83, just outside our redshift window. Given typical photo-z uncer- 
tainties, a significant fraction of objects from this cluster still leak into our 
redshift selection window (see Fig. l 111 . 



Vmax is the total volume in which the galaxy could be located and 
yield an [O II] flux above our flux limit at the observed wavelength 
and yield a A'-band magnitude between our bright and faint K- 
band limits. 



We calculate Vmax for the i th galaxy as 



V m 



= n 



dz 



(12) 



where Q is the angular area of our survey (107.8 arcmin 2 
for CDFS and 29.1arcmin 2 for FIRES) , dV c /dz is the differen- 
tial comoving volume (e.g., iHogjl 19991) . and F viSt i(z) is a visi- 
bility function for each galaxy which is 1 when both its AT-band 
and [O II] fluxes are above our survey limit at that redshift (wave- 
length), and otherwise. The calculation of the two terms (based 
on the AT-band and [O II] flux limits) making up F v i 3t i(z) are out- 
lined below. For reference, a galaxy which is visible at all redshifts 
(from z = 0.889 to zi = 1.149) in ROLES would be drawn from 
a volume of 7.4 x 10 4 Mpc 3 . 



3.4.1 K-band flux limit 

For a galaxy observed at redshift z obs with if-band magnitude 
K obs , the JC-band magnitude at a given redshift z is estimated to 
be: 



K z = K obs + 5log(d oba /d z ) + (k-c 



k-corr obs ) (13) 



where d obs and d z are the luminosity distance at the observed red- 
shift and the redshift at which we wish to estimate K respectively. 
k-corr obs and k-corr z are the fc-corrections at the same redshifts. 
For simplicity we adop t the if-band differential fc-correction of 
lGlazebrooketal.Nl995h . which is based on an average model SED 
for normal galaxies: 



k-corr(z) — 



-2.58z + 6.67z 2 - 5.73z 3 - 0.42z 4 



(14) 



2.362 + 3.822 2 - 3.532 3 + 3.35z 4 

In any case, in the A"-band the differential fc-correction de- 
pends only slightly on the galaxy type (for z<1.5), and due to our 
narrow redshift slice the fc-correction is constant to a good approxi- 
mation. In practice, the effect of applying a fc-correction or not only 
makes a ~1% difference to the total volume probed. 



3.4.2 [Oil] flux limit 

Our [O II] line flux limit is not simply a constant flux limit (as is 
sometimes assumed in similar surveys), but is in fact a function 
of wavelength due to the numerous bright night sky lines. Since 
we have propagated an estimate of the noise in our spectra, Nij, it 
is a simple matter to calculate the average flux limit as a function 
of wavelength. We calculate the average noise spectrum for each 
mask and then take the shallowest mask as a conservative limit. 
The actual value we use is 4Nij, since we adopt a 4a detection 
threshold. This refers to the flux in the unsmoothed image, which 
is what is used to make flux measurements. Since object detection 
is performed in the gaussian-smoothed images, some detections can 
appear below the formal flux limit in the unsmoothed frame (as it is 
only by smoothing to increase the S/N that they pass the detection 
threshold). Fig. Q3] plots the fluxes of our [Oil] lines versus their 
observed wavelengths. The solid line shows our 4a [O II] limit. As 
can be seen, the limit varies from approximately 3 x 10 -18 erg s _1 

17 —1 —2 

erg s cm 



cm 2 in the deepest regions to around 1 x 10 
underneath the brighter sky lines. Those few points which appear 
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Figure 13. Flux versus wavelength for the lines with Pou > 0.5 (i.e. 
most-likely [O II] detections, solid circles) and Poll < 0.5 (open circles). 
The solid lines show the 4cr flux limit derived from 4 X the flux measure- 
ment error determined from the average noise spectrum for each mask. The 
limit used in ROLES is that of the shallowest mask (upper line) to be con- 
servative. Circles show CDFS detections and squares show those from the 
FIRES field. Points may appear below the detection threshold, since the flux 
measurements are made from the unsmoothed data frames; whereas the de- 
tection limits refer to the noise level from the Gaussian smoothed frames. 
The few points which are above the detection limit but below the measure- 
ment limit are rejected from further analysis. 



below the flux limit of the shallowest mask (upper solid line) are 
rejected from further analysis. 

Averaged over all of our [O II] lines, calculating the total 
Vmax in this way using the detailed noise spectrum gives a value 
20% lower than had we calculated it assuming a constant flux limit 
of the flux in our faintest line. 



3.5 Higher mass comparison samples 

In order to extend the mass range considered in this work to higher 
masses, we employ two external samples. The first is the Gem- 
ini Dee p Deep Survey ( GDD S). Details of the GDDS data are 
given in Abra ham et all J2004I) . Briefly, this is a if-selected sur- 
vey at somewhat brighter magnitudes (higher stellar masses) than 
ROLES, but still much fainter than typically attempted by spectro- 
scopic surveys. GDDS used N&S with the GMOS spectrograph on 
the 8-m Gemini telescope. Spectroscopic exposure times were typ- 
ically around 20 hours. Four fields totalling an area of 121 arcmin 2 
were observed, but the regions selected for spectroscopy were nor- 
malised to a larger photometric survey area of 554.7 arcmin and 
so this is the area that should be used when considering sampling 
volume, cosmic variance, etc. For the GDDS fields, the photome- 
try available is VRIzK, with 2/4 fields lacking i?-band. They use 
colour selection (/ — K) to preferentially target non-star forming 
galaxies, but include a sparse (~1 in 6) sampling of bluer (star- 
forming) galaxies. Redshifts are determined interactively, by visu- 
ally comparing spectra to templates. This likely yields accurate red- 
shifts where multiple features are visible, but may, for example, not 
be as highly complete for low SFRs as ROLES (for example, a sin- 
gle 4cr emission line in GDDS is less likely to be spotted by eye as 
in our automated search for this specific class of object). 

For the GDDS data, the photometry was fitted to a similar 
grid of models as for deriving the ROLES stellar masses and SFRs 



(ED- One minor difference i s the dust extinction law used in the 
model fits (SMC for GDDS vs lCalzetti etafll2000l for ROLES). 

The second comparison sample comes from the ES Q public 
FORS2 spectroscopy in GOODS-S JVanzella et al]|2008l) . This is 
largely a z-band selected survey, but with some colour-selected 
(including photometric redshift-selected) subsamples. This survey 
uses 1" slits and is flux calibrated based on the broadband photome- 
try of continuum sources. Since this sample covers exactly the same 
region of sky as ROLES (and thus the same multiwavelength data 
are available), it provides the ideal sample to extend our survey to 
higher masses. Our flux calibration is compared with spectroscopy 
in common with this survey in i]2.5.1l Exposure times with FORS2 
(on the 8-m VLT) were typically around 4 hours, and so the spec- 
troscopic data are comparable in limiting [O II] flux, or deeper than 
the ROLES data. We take spectroscopic redshifts, confidence flags 
and ID flux calibrated spectr J*^l for galaxies which directly overlap 
with our observed LDSS3 pointings in CDFS. Only the two highest 
confidence classes of redshifts (A and B) and objects in the ROLES 
redshift window (0.889 < z ^ 1.149) are used. From the ID spec- 
tra, [O II] luminosities are measured in a way analogous to the 
ROLES method. From measurements of many faint emission lines, 
the 4.0cr flux limit is estimated to be 6 x W~ is erg s' 1 cm~ 2 . Ob- 
jects are matched to our photometry and stellar masses fitted, com- 
pleteness calculated and [O II] SFR estimated in exactly the same 
way as for ROLES. The (targeting) completeness for the FORS2 
sample is around 40% in each A"-band magnitude bin, except for 
the brightest bin (highest mass galaxies) which is around 70% com- 
plete. The success ra te of measuring redshifts is taken to be 72% 
dVanzella et alj|2008l) independent of magnitude (which is likely a 
reasonable approximation for emission line galaxies). 



4 RESULTS & DISCUSSION 
4.1 [O II] luminosity function 

Fig. [14] shows the [O II] luminosity function (LF) measured by 
ROLES. For the i th bin of log luminosity, the number density, $i, 
is calculated using the \IVmax method as 

PoiX.i 



^ V m 



(15) 



where Pou is the weighting due to the probability that the emis- 
sion line is [O II], V ma x,i is the maximum volume in which galaxy 
i could be located and have been found in the ROLES survey, 
and Wi is the weighting given by the spectroscopic completeness 
(eqn.ll It. Error bars are calculated by s caling Poisson distributions 
following the method o fZhu et all J2008l) . as summarised below. A 
scale factor within each log luminosity bin is determined from an 
effective weight, W e ff, 



W ef f = 



( Vmax ) i 



and the effective number of objects, N e ff, 



Neff = 



T — 

^ (V™. 



a x ) 1 



W. 



eff- 



(16) 



(17) 



15 see http://www.eso.org/science/goods/ 
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Figure 14. [O II] luminosity function. Filled circles and open squares show measurements made in the CDFS and FIRES fields respectively. Error bars are 
derived from Poisson errors only. In the left panel, open triangles and diamonds show the measurements o f lzhu et alj 1200 8|) in their 0.752 < z < 0.926 and 
0.926 < z < 1.099 bins, respectively. The dashed lines show the best fit line found bv lZhu et alj 1200 8|) for log £[o/j] > 42 and their assumed faint end 
slope of a = —1.3 at low luminosities, renormalised to pass through our data. The thick line shows our best fit faint end slope of a = —1.5. Red open circles 
and blue asterisks show comparison data from the higher mass samples of GDDS and FORS2 (offset slightly for clarity) respectively. In the right panel, the 
z~l ROLES measurements are now compared with the z~0. 1 measurements made in SDSS (for mass cuts similar to the ROLES' mass range: diamonds; and 
for no mass cuts: thicker diamonds with error bars). 



The upper and lower limits are then ca lculated for th e effective 
number assuming a Poisson distribution dGehreld[l986l provides a 
useful way to approximate this) and these numbers are scaled by 
W e f f to give the upper and lower SFRD limit in each mass bin. 

Fig. [14] shows the measurements for the CDFS and FIRES 
fields separately (filled circles and open squares respectively). 
Since cosmic variance is likely to be important, and is not included 
in the error bars, the difference between the two fields gives an 
indication of the magnitude of thi s effect. Estimating the square 
root of cosmic variance, av, from lSomerville et all ( 2004) using 
the number density of ROLES galaxies in these fields gives a value 
of o y ~ 0.7. This is consistent with the observed difference be- 
tween the two fields. 

Ope n triangles and di amonds (left panel) show the measure- 
ments of IZhu et alj j2008l) in their 0.752 < z < 0.926 and 
0.926 < z < 1.099 bins, respectively (error bars omitted for clar- 
ity). They suggest that the z~ 1 [O II] LF is best fitted by a dou- 
ble power law model with a turnover around log(.Lrojj]) ~ 42. 
Their sample was insufficiently complete at lower luminosities to 
constrain the faint end slope, a (dN/dL oc L a ), so they assumed 
a = —1.3. Their fit and assumed slope are shown as the dashed 
lines in Fig. [14] The faint end line has been normalised to fit 
through the ROL ES data. It can b e seen that the ROLES data do not 
overlap with the lZhu et alj J2008) sample brighter than the turnover. 
This is not unexpected, since ROLES is designed to be a low stellar 
mass sample, and thus low mass galaxies cannot su stain the high 
SFRs implied by such high [O II] luminosities (e.g., Noes ke et al.1 
120071) . This mass selection is likely to be the biggest difference in 
comparing the LF between the two samples, and may be responsi- 
ble for the apparent offset in normalisation where the two samples 
overlap in luminosity. We fit a power law to the ROLES data and 
find a faint end slope of a ~ —1.5, although this is not well con- 
strained. The main purpose of performing this fit is to obtain an 
indication of the turnover due to incompleteness in [O II] luminos- 
ity. This appears to occur around log(Lrojn) ~ 40.5, assuming 
the intrinsic distribution can be well approximated by a power law. 
Using eqn. [9] it can be seen that SFR of IMq yr~ l corresponds 



to a log luminosity of 40.7, assuming 1 magnitude of extinction at 
Ha. This suggests that ROLES is becoming incomplete in SFR at 
SFR~ O.7M0 yr^ 1 . However, this amount of extinction is likely 
an overestimate for these low mass galaxies (see e.g., G10, and next 
section), and thus the limiting SFR is likely lower. 

In the right panel, the ROLES data are now compared with 
the z~0.1 [Oil] LF measured in the SDSS (G10). Open dia- 
monds (error bars omitted for clarity) show the LF for a subsam- 
ple cut in stellar mass to approximately match the ROLES range 
(8.5 < log(M/M0) < 9.5), and thicker diamonds show the 
whole sample of SDSS galaxies, unrestricted in stellar mass. The 
mass selection obviously makes a large difference to the shape of 
the [O II] LF, as mentioned above. However, if we compare the 
ROLES data to the SDSS sample unrestricted in mass, we find 
that the measurements (where the samples overlap, at the faint end) 
agree within the errors. Assuming at z~l ROLES' mass galaxies 
dominate the [O II] LF in this luminosity range, this is a reasonable 
comparison. The fitted faint end slope in the SDSS data is a faint = 
-1.67 ± 0.02, compatible with the ROLES fit. At the bright end 
in SDSS, airight = —2.83 ± 0.29, somewhat shallower (but com- 
patible wi t hin th e uncertainties) than the z~l measurements from 
IZhu et alj J2008t) which are around -3.2 - -3.0 ±0.1 for the data 
shown in Fig. [14] Although the bright and faint end slopes are con- 
sistent with their z~l values, the position of the turnover clearly 
evolves. The measured values are log Lto = 4 1 .0 in SDSS (G10) 
and log Lto = 41.9 at z~l jzhu et al.ll2008l) . IZhu et al.1 d2008t) 
propose a parameterisation for the evolution of the bright end of 
the [O II] LF (which they measure between 0.75 < z < 1.45) by 
considering the the redshift evolution of objects where the number 
density is 10~ 3 5 dex' 1 Mpc~ 3 as logL [0/I] = 0.46z + 41.85. 
Extrapolating this relation to z~0. 1 would give an expected lumi- 
nosity for this number density of objects of log Lyon] = 41.9. In 
fact, the value measured from SDSS is closer to log Lyon] = 41.4, 
half a magnitude brighter. Thus the evolution is stronger between 
z~0.1 and z~l than that expected from a simple extrapolation of 
that between 0.7 < z < 1.5. 
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Under our simple conversion of [O II] luminosity to SFR, 
Fig. [14] can be viewed as a SFR function. This may suggest that 
the SFR function has remained constant in shape but has simply 
shifted toward lower SFR (lower [O II] luminosity) at lower red- 
shifts. This could imply that galaxies have simply decreased their 
SFR by the same amount, on average, over this redshift range. The 
[O II] LF is a relatively blunt tool for understanding how the SFR 
of the general population evolves, and so in the next section we ex- 
amine the star formation rate density as a function of redshift and 
stellar mass. 



4.2 Star formation rate density 

For each mass bin, the star formation rate density, psfr, is calcu- 
lated using the l/Vmax method as 



PSFR 



(18) 



where the symbols have the same meaning as for eqn.[T5]and SFRi 
is the star formation rate (derived from Lyon] following eqn.[5]l of 
the i th galaxy. Error bars are again calculated by scaling Poisson 
errors, as for the [O II] LF, but this time in mass bins instead of 
luminosity bins. 

Again, the results are plotted separately for the two fields to es- 
timate the effect of cosmic variance. Slightly different-sized mass 
bins are used in the two fields due to the different sample sizes. 
Bin sizes are chosen to give comparable numbers of galaxies in 
each mass bin. SFRDs, errors, and the numbers of objects used 
are tabulated in Table [2] We note that the mass cuts imposed re- 
move a few objects from the initial sample summarised in §2.5 
(and used for the LF analysis). For objects passing the flux limit 
and magnitude cuts, we have 212 in CDFS (17 of these are given 
fractional weights, i.e., have 0. 1 < Poii <0.9) and 90 in FIRES (40 
of these with fractional weights). Once the mass cuts are imposed, 
these numbers are reduced to 199 for CDFS and 86 for FIRES, 
split by mass as tabulated in Table[2] For a higher mass comparison 
sample, we also include data from GDDS and FORS2. We have 
recalculated SFRDs for GDDS using their raw SFRs, complete- 
nesses and weightings (kindly provided by S. Juneau) in exactly 
the same way as for ROLENS Results for our default [O II] SFR 
estimator (Eqn.[9} are shown in the top left panel of Fig. [15] This 
is the same transformation from [O II] luminosity to SFR as used 
in GDDS (J05). Filled circles and open squares show the results 
for the CDFS and FIRES fields respectively. The slightly differ- 
ent mass ranges probed by the CDFS and FIRES fields are due 
to the different intrinsic colour distributions of the fields and the 
fixed A"-band limits used, as described in i)3.1| Horizontal error 
bars show the range of the mass bins, with the data points indi- 
cating the centre of the adopted bin. Asterisks indicate the median 
mass within the bin. GDDS points are shown as open red circles. 
The fact that the GDDS and FORS2 points seem to join smoothly 
onto the ROLES points at higher masses is a useful independent 



16 For GDDS, we use data in the redshift range 0.88<z^l.40. The mass 
bins are chosen to give comparable numbers of objects in the overlap with 
ROLES' redshift range (0.88<z^l.l5) and using data in just this range 
makes negligible difference to our [O II]-based SFRD results. However, 
UV SFR estimates are only available at z>1.2 (due to their limited very 
blue photometry) and so we extend the redshift range slightly higher than 
this to allow comparison of UV SFRDs, in the next section. 



Table 2. Nominal [O II] SFRD in bins of stellar mass. M h M u and M m 
refer to the lower and upper limits of the mass bin (lcr Poisson error) and 
the median log(mass) within the bin, respectively. SFRD, pg pn, is given in 
units of 10~ 3 A/q yr~ 1 Mpc~ i per dex of stellar mass. Again, subscripts 
u and I denote upper and lower limits. N ga i gives the number of galaxies 
in the mass bin (for GDDS, numbers in parentheses refer to the subsample 
with UV SFRs: 1.20 < z < 1.40 only). 
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check of our absolute flux calibration and means that our calibra- 
tion is likely as good, or better, than our estimated 30% absolute un- 
certainty. For reference, the two diagonal crosses show the results 
(error bars omitted for clarity) from the smaller subset of ROLES 
data presented in paper 1, where we first observed the turnover in 
the SFRD towards these lower mass galaxies. It can be seen that 
our much larger dataset now allows a better determination of the 
shape of the SFRD, and the use of two separate fields allows the 
effect of cosmic variance to be assessed. 

The solid curve shows the same measurement made at z~0.1 
in the SDSS (G10). The dashed curve is the z~0.1 SFRD renor- 
malised by a factor of 3.5 to match the amplitude of the z~l data. 
It can be seen that this simple renormalisation leads to good agree- 
ment between the shape of the z~l and z~0.1 SFRD, suggesting 
that simple density evolution may be enough to explain the differ- 
ence in SFRD between the two epochs. 



4.2.1 Limits on the SFRD from ROLES galaxies with no 
[O II] detection 

So far, only galaxies with emission line detections have been con- 
sidered in the analysis. However, galaxies may be present in the 
ROLES redshift (and mass) window, but since we do not detect 
an emission line, we can assign neither a redshift nor a SFR. In 
order to estimate how much these objects may contribute to the 
z~l SFRD, we perform the following test. We repeat our analy- 
sis for all remaining galaxies with no emission line detections. We 
exclude objects known to lie outside our redshift window from se- 
cure redshift measurements (from public spectroscopy or our own 
multiple line detections), but use the photo-z's of all remaining ob- 
jects to assign a probability, Poii, that they lie in the ROLES red- 
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Figure 15. The SFRD as a function of stellar mass using various SFR estimators. In each panel the ROLES data are shown separately for the CDFS and FIRES 
fields as filled circles and open squares respectively. (Diagonal crosses show the mean SFRD in the subsample of CDFS from paper 1, for reference.) The 
higher mass GDDS points are shown as red open circles and FORS2 as blue diamonds. The equivalent SFR estimator at z^O.l taken from SDSS is shown as 
a solid curve. The dashed curve is the solid curve renormalised by a factor of 3.5 to agree with the observed z~l [O II] SFRD. The upper left panel shows 
our default [ O II] indicator (constant extinction); upper right panel shows [O II] SFR corrected by observed B-band luminosity following the prescription of 
Moustakas et al. 12006); lower left panel shows the empirical mass-dependent correction of G10 to [O II ] and the lower right pan el shows the SFRD obtained 
by SED-fitting. Orange open triangles show the UV-determined SFRD for the spectroscopic sample of[Cowie & Barger (2008). In the lower left panel, the 
open diamonds show [OII]+24/im SFRs from DEEP2 and the dotted line shows the B aimer-decrement corrected Ha SFRD (also scaled by a factor of 3.5) 
from SDSS. See text for discussion. 



shift window. For simplicity, we only do this for the CDFS field. 
The much greater area and spectroscopic follow-up means that this 
provides better statistics than the FIRES field, anyway. This yields 
232 additional galaxies with non-zero probabilities of being in the 
ROLES redshift range. For these objects, we take the peak of their 
photo-z PDF as their true redshift, and measure the 4<r upper limit 
on the flux at the expected position of [O II] at this redshift. Ob- 
viously, due to the size of the typical photo- z error (Az ~0. 1, 
iForster Schreiber et al.l2o6o1IWuvts et all2 008 ) , this location is not 
an accurate estimate of the position of the expected [O II] line, but 
provides a representative estimate of the noise and thus an approxi- 
mate 4a upper limit to the [O II] luminosity. Repeating this test us- 
ing random sampling of the photo-z PDF instead of the peak yields 
similar results. These upper limits are then added to the CDFS cat- 
alogue and propagated through the SFRD analysis. The final anal- 
ysis including these limits results in an increased SFRD in each 
mass bin which approximately coincides with the upper la pois- 
son error bar (as plotted in Fig.ll5t in all cases. These limits are not 



shown on the plot to preserve the clarity of the results from the de- 
tected galaxies. Thus, ROLES does not miss significant star forma- 
tion from galaxies which we do not detect down to our [O II] flux 
limit. 

Additionally, we should also note that we have not explicitly 
removed AGN from our sample. The exclusion of AGN will only 
lower the [O II] SFRD in ROLES. Examination of the MIR colours 
of ROLES galaxies and X-ray point sources using the extensive 
multiwavelength data in CDFS shows that AGN contamination is 
of extremely minor importance in our survey. 

4.2.2 Empirically correcting the [O II] SFR 

So far, a simple constant scaling between [O II] luminosity and 
SFR has been assumed. As discussed in G10, this is too naive and 
overestimates the SFRD at low masses and underestimates it at high 
masses. Assuming the mass dependence of [O II] luminosity on 
SFR (from metallicity, dust extinction, ionisation parameter, etc.) is 
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the same at z~l as z~0.1, any trends with redshift are still robust, 
but the overall shape of the SFRD as a function of mass is likely 
ske wed by this simple estim ator. 

iMoustakas et all ( 2006) present an empirical method for cor- 
recting [O II]-SF Rs using rest-fr ame B-band luminosity (which is 
often used, e.g jzhu etal]|2008l) . Although G10 showed that this 
method underestimates the size of the correction needed at high 
stellar masses, we consider it for completeness and to allow com- 
parison with other works using this procedure. Synthesizing rest- 
frame B-band absolute magnitudes, Mb, from the SED-fitting pro- 
cedure used to derive stellar masses ( ij3.lt we apply this correction 
to our [O II]-SFRs for each galaxy and the resulting SFRD is shown 
in the top right panel of Fig. Q3] It can be seen that this correc- 
tion decreases the SFRD in low stellar mass objects and increases 
the SFRD in the highest stellar mass objects. The solid line shows 
the identical measurement made at z~0.1 in SDSS (G10) and the 
dashed line shows this line renormalised by the same factor (3.5) 
which gave agreement for the constant luminosity [O II] SFRD. 
Using this alternate estimator, all the ROLES points still give ap- 
proximately the same agreement with the scaled local SFRD as for 
the default [O II] estimator, however only one of the three GDDS 
points and one of the four FORS2 points are now consistent with 
this curve at the lcr level. 

In the lower left panel of Fig.[T5] we plot the SFRD calculated 
using [O II] with an empirical mass-dependent correction derived 
by G10 to reconcile [O II] SFRs with Balmer decrement-corrected 
Ha SFRs at z~0.1. This assumes that the mass-dependence of 
[O II] luminosity scales in the same way at this redshift as locally. 
Again solid and dashed curves show the corresponding SDSS mea- 
surement and the same, renormalised by a factor of 3.5, respec- 
tively. A fairer way to make the comparison with the local sam- 
ple is using the Balmer decrement-corrected Ha SFRD rather than 
the [Oil] SFRD since, as shown in G10, the [Oil] flux limit in 
SDSS is not low enough for the [O II]-based SFRD to have con- 
verged. Using the SFR functions in mass bins, G10 showed that 
ROLES reaches low enough [O II] fluxes for the SFRD to have 
converged (assuming the shapes of the z~l SFR functions in bins 
of mass follow the same shapes as at z~0.1). The dotted line shows 
the z~0.1 Balmer decrement-corrected SFRD scaled by 3.5. This 
scaled Ha curve passes through the z~ 1 data systematically higher 
than the [O II] curve, for the reason just mentioned: incompleteness 
in the z~0. 1 [O II]-selected SFRD, which becomes important at 
higher masses. For this reason, we choose to renormalise the scale 
factor between the z~l and z~0. 1 (Ha) data. A factor of 2.6 is a 
better fit than the 3.5 previously found from the (incomplete z~0. 1) 
[O II] data. From now on, we will use this factor when comparing 
the empirically-corrected [Oil] SFRD with the z~0.1 Ha SFRD. 
This curve will be plotted in the next figure. 

The fact that the renormalised local measurement agrees with 
our z~l measurement is largely a restatement of the agreement 
for our nominal [O II] measurement (upper left panel). However, 
since this estimator is our preferred estimator (using [O II]) of 
total SFR^l we can now compare with z^l "t otal" SFRs using 
[OII]+24/im SFRs from IConselice et all ( l200l DEEP2), shown 
as red open diamonds. Although the error bars are large and the 
overlap with the GDDS, FORS2+ROLES sample in mass is min- 
imal, the agreement at log(M*/Mo)~10 is encouraging. Now 



17 The SFRD data points may straightforwardly be obtained by applying 
the G10 correction (SFR corr = SFRo/{-lA24tanh[(logM/M Q )- 
9.827)/0.572] + 1.700}) to the (SFR ) values in Tableff] 



comparing with the rescaled z~0. 1 data, we see that the DEEP2 
points are both consistent with the scaled local value. 



4.2.3 SED-fit SFRs 

In the lower right panel of Fig.[l5] we examine the SFRD where the 
SFR of each galaxy is estimated from SED-fitting of the deep, mul- 
tiwavelength photometry at each galaxy's spectroscopic redshift (as 
used for the stellar mass-fitting in 5j3j}. This method is now com- 
pletely independent of [O II] luminosity (except that for ROLES, 
a galaxy must have a minimum [O II] luminosity to be selected in 
the survey). This is largely equivalent to a rest-frame UV luminos- 
ity SFR (with the longer wavelength information constraining the 
contribution of dust extinction). Hence, for the local SDSS com- 
parison sample, we use the u-band SFRD with the dust correction 
estimated from the Balmer decrement (G10). Again, the solid line 
shows this, and the dashed line is the same renormalised by 3.5 
in number density. At low stellar masses (ROLES), the "UV" and 
[O II] SFRDs approximately agree (with the scaled local predic- 
tion and hence with each other), albeit with larger scatter compared 
with [Oil], and some suggestion of a systematic offset, since all 
but one of the SED-fit points lie above the line. The larger scatter 
is likely due to the significant photometric errors on such faint ob- 
jects. It should be noted that in order to obtain UV SFR estimates 
for GDDS, the redshift range is restricted to z> 1 .2. In order to min- 
imise the effect of evolution (which becomes significant between 
z~1.0 and z~1.6, see fig. 2 of J05), we restrict the redshift range 
to 1.2 < z 1.4. This results in very small samples of galaxies in 
the two lowest GDDS mass bins (2 and 6) and a modest sample of 
16 in the highest mass bin. No such extension of the redshift range 
is necessary for the FORS2 data, and this sample covers exactly the 
same redshift range as ROLES, with many more object in FORS2 
possessing UV photometry than in GDDS (20-30 vs 2-6). 

For comparison, open orange triangles show estimates of 
the extinction-corrected UV SFR from the spect roscopic redshift 
(0.9< z ^1.5) survey of ICowie & Bargeil 12008), transformed to 
our IMF. These UV SFRD estimates agree well with the SED- 
fit SFRD from the FORS2 data. For the log(M*/M Q )~ll mass 
bin where the GDDS sample size is reasonable, the agreement 
is also good. The lower mass GDDS bins measure significantly 
lower SFRDs than ICowie & Barged J2008f). Some ca ution is re- 
quired when comparing to the lCowie & Barged ( 120080 data, since 
their UV SFR is not measured in exactly the same way as in 
ROLES or GDDS. ICowie & Barged 12008) derive a slightly dif- 
ferent calibration of the UV SFR than that used here. They also use 
NIR luminosity to estimate the contribution of older stellar pop- 
ulations to the UV l uminosity and subtract this before converting 
to SFR. Given that ICowie & Bargeil (2008) derive these correc- 
tions to the UV SFR self-consistently within their data, we do not 
attempt to correct for the slightly different overall normalisation 
they use in UV luminosity-SFR, and we just apply a correction for 
the change from their Salpeter to our BG03 IMF. Undoing b oth of 
these effects would only increase the ICowie & Barged d2008h UV 
SFRD, strengthening the excess seen at higher masses relative to 
the GDDS data. In light of the fact that the sample size in these bins 
is tiny for GDDS, and that the FORS2 data (base d on much larger 
numb ers of objects) give better agreement with the lCowie & Barged 
120081 results, we feel justified in rejecting the GDDS SED-fit SFRD 
values from further analysis. 
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Figure 16. Best estimates of SFRD. Symbols indicate the different surveys, 
as for Fig. 1151 Colours refer to the method of estimating SFR: black symbols 
indicate [O II]-based measurem ents, blue symbols SED-fit SFRs, orange 
symbols ICowie &~B arger 1 200§)'s UV measurements, and red points are 
for the DEEP2 [O II]+24/im data. The solid line shows the z~0. 1 Balmer- 
decrement corrected Ha SFRD, the dashed line is this scaled by 2.6, and the 
dotted line this scaled by 6. The UV/SED-fit based SFR indictors predict a 
slightly higher overall normalisation for the z~l SFRD than the [O II] es- 
timates, but the overall shapes are the same. See text for discussion. 



4.2.4 Best estimates of the z~i SFRD 

After discussing the various different SFR indicators which may 
be used with these datasets, it is useful to collate all the results of 
the best estimators. Our preferred estimators are the empirically- 
corrected [Oil] estimator of G10, and the SED-fit SFR estimates 
(after dropping the GDDS results). These are compiled into one 
plot in Fig. [16] Symbols indicate the different surveys, as for 
Fig. [15] The colours refer to the method of estimating SFR: black 
symbols indicate [O II]-b ased measurements, blu e symbols SED- 
fit SFRs, orange symbols ICowie & Barged |2008)'s UV measure- 
ments, and red points are for the DEEP2 [O II]+24/im data. The 
solid line shows the z~0.1 Balmer-decrement corrected Ha SFRD 
and the dashed line is this scaled by 2.6. Viewing all this data 
together, colour-coded by SFR estimator, it is now clear that al- 
though the [O II]-based SFRDs seem consistent with a renormali- 
sation of ~2.6 times higher than the local SFRD, the SED-fit and 
UV-based estimates are systematically somewhat higher. The dot- 
ted line shows the local SFRD renormalised by a factor of 6, which 
is closer to the z~l value for the UV/SED results. 

Both datasets ([O II] and UV) show comparable shapes for 
the mass-dependent SFRD - shapes similar to that determined at 
z~0. 1 from Ha measurements - but the different estimators dis- 
agree somewhat in the overall normalisation for this curve. We 
briefly discuss reasons for this possible disagreement (although it 
should be noted that uncertainties between SFRDs measured with 
different indicators at these redshifts are typically as large or larger 
than we see here, e.g., Hopki ns & Beacom 2006). 

UV estimate methodology - G10 compared UV (u-band and 
FUV)-based estimates of the local SFRD with those of [O II] and 
Ha and found that they gave comparable estimates (certainly not 
discrepant systematically in normalisation by the factor of ~ 1.7 
seen here). Although the method of SED-fititng used here differs 
somewhat fr om the dust-corre c ted res t-frame UV luminosities used 
in G10, the|c owte &Barger (2008) data use a very similar ap- 



proach to estimating UV SFRs as that used in G10. So, the fact 
that the Cowie & Barget] d2008l) and SED-fit methods give compa- 
rable SFRDs suggests that this fitting approach is not responsible 
for the differences between [O II] and UV SFRDs at z~l. 

Aperture corrections - Both the photometric (UV/SED-fit 
SFR) and emission line ([OII]-SFR) estimates require aperture 
corrections, as discussed in ij2.5.1l i j3.ll respectively. The former 
use the A'-band light scaled from the smaller colour aperture to 
the value u sed for the total mag nitude; and the latter use the cal- 
ibration of Vanzella et al. (2008), assuming our z~l galaxies are 
seeing-dominated in our spectroscopy. Again, we can use the fact 
that the data measured using our method (e.g., the FORS2 data us- 
ing our photometry an d SED-fitting technique compared with the 
ICowie&Bargej|2008l data) agree well with other results in the 
same mass range as a check that our aperture corrections cannot 
be egregious. Similarly, the agreement between the GDDS [O II]- 
SFRD and our [O II]-SFRD (for the highest mass galaxies where 
the aperture corrections should be the largest) shows that uncer- 
tainties in the aperture corrections are relatively unimportant, since 
GDDS do not explicitly use an aperture correction and for the 
FORS2 data we used the correction based on ROLES data. 

Incompleteness in [O II] SFRD - It is worth emphasizing that 
the lower overall SFRD measured from [O II] versus UV is un- 
likely to be due to missing flux below the [O II] limits of the var- 
ious surveys. As we have estimated in i)4.2.1| galaxies in ROLES 
with no detectable [O II] emission, but plausibly within our red- 
shift window would, at most, move the [O II] SFRD to the top of 
our 1<t SFRD errors, causing better agreement between [O II] and 
UV SFRDs. However, two important points must be borne in mind. 
Firstly, this only applies for the ROLES dataset, which is [O II]- 
selected. For FORS2 and GDDS, no such explicit dependence on 
[O II]-selection applies. Secondly, since ROLES is [O II]-selected, 
including more objects in the [O IIJ-SFRD would also increase the 
contribution to the UV SFRD, since these new entries would then 
have their SEDs fitted in the same way as for the galaxies currently 
in ROLES, driving up the UV SFRD. Additionally, the checks done 
in G10 show that the SFRD has likely already converged above the 
ROLES' (and other surveys employed here) flux limit. 

IMF - The [O II] and UV luminosity trace SFR via their de- 
pendence on the presence of high mass stars. However, they are 
sensitive to stars in dif ferent mass rang es: ;> lOM© for [O II] and 
> 5M© for UV (e.g., Kennicutt 1998). Thus, a different slope of 
the stellar IMF can lead to a different normalisation between the 
two estimators. The agreement between the estimators at z~0.1 
would mean that, if the IMF is responsible, then it would have to be 
an evolutionary effect (which affects galaxies of all masses approx- 
imately equally). An evolving IMF has been suggested as a way 
to reconcile estimates of the growt h of stellar m ass with the cos- 
mic star formation rate history (e.g., Dave 2008 - an evolving high 
mass turnover which moves towards lower masses at higher red- 
shift). Another possibility is that instead of a universal IMF which 
evolves, the IMF is different (but constant with time) for starbursts 
versus more quiescent star formation. There is good motivation for 
assum ing that the IMF in starbursts may be relatively top heavy 
(e.g.. lLaceT et al. 2008 and references therein). Then, the evolving 
contribution of starbursts to the total SFRD would lead to an evo- 
lution of the relative normalisation of the indicators. It should be 
noted that, if the incidence of starbursts is higher at higher redshift, 
then an increasing contribution from a top-heavy IMF at higher 
redshift would go in the wrong direction, i.e., one would expect 
that the indicator more sensitive to the most massive stars ([O II]) 
would show a higher SFRD than the other (UV), under the assump- 
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tion of the incorrect IMF. Conversely, the Dave ( 2008) model for 
IMF evolution goes in the correct direction to explain our data. 

Evolution of the empirical [O II] correction - We have cor- 
rected our SFR based on [O II] luminosity as a function of mass 
using the local empirical calibration of G10. This accounts for 
the (mass-dependent) trends of [O II] luminosity-SFR on metal- 
licity, dust and ionisation parameter. If the z-~l trends with mass 
are still comparable for all of these parameters, then the relation 
is still valid. However, it is expected that parameters such as dust 
and average metallicity will evolve wit h cosmic time. For ex ample, 
for the higher mass galaxies in GDDS, ISavaglio et ail fe005h found 
that the metallicity of log(M*/M Q )~9.5 galaxies at z~0.7 is ~0.1 
dex lower than local galaxies of the same mass. If metallicity evo- 
lution were responsib le, differential evolution as a function of mass 
( Savaglio e t alj|2003) should lead to a differential correction as a 
function of mass, which would change the shape of the SFRD. Of 
course it is possible that some conspiracy between dust and metal- 
licity evolution changes the nor malisation but not th e shape. The 
metallicity dependence of [O II] jKewlev et alj|2004h is such that 
the predicted correction to the [O II]-SFR would be lower for a 
decreased metallicity, i.e., that our [O II] SFRD should be overes- 
timated relative to the UV, if this is the cause. In practice, we have 
corrected empirically for not just the metallicity, but also the dust 
extinction, etc. Assessing the validity of this empirical correction at 
z~ 1 will have to await further follow-up data, such as Ha/H/3 spec- 
troscopy of ROLES galaxies, and such work is underway. Even if 
the empirical correction is not correct in detail at z~ 1 , applying this 
correction is still better than not applying any mass-dependent cor- 
rection. The fact that the shape of the [O II]-SFRD is similar to that 
of the UV SFRD, suggests that evolution in the overall correction 
might not be responsible. 

We have attempted to carefully unify measurements of the 
SFRD from different surveys (which probe different mass ranges) 
at the same epoch (z~l), in an attempt to study trends with stel- 
lar mass. Similar caution must be taken when comparing different 
surveys at different redshifts in order to study evolution. It is worth 
noting that we now disagree with the result of J05 who found that 
their intermediate mass bin (10.2 < log(M*/M )< 10.8) reached 
the local value of the SFRD by z~l. Our result using their data 
shows that the SFRD for this mass of galaxies is at least 2.6 times 
the l ocal value. This differenc e is because J05 used the Ha data 
from lBrinchmann et al.1 d2004l) to estimate the z~0 SFRD and thus 
has a very different shape from the [O II] SFRE{3 (G10), scaled 
to their IMF. As can be seen from comparing the top left panel of 
Fig. [T5]( which is the [O II] estimator J05 use) with the lower left 
panel (see also G10), the local SFRD (solid curves) derived using 
the different estimators is significantly different in the GDDS mass 
range. 



5 CONCLUSIONS 

We describe the methodology for ROLES, a survey for dwarf galax- 
ies at z~ 1 . Using [O II] A3727 we have estimated the star formation 
rates (SFRs) of galaxies with spectroscopic redshifts, and masses 

18 Also compare the shapes of the dotted curve in the lower left panel of 
Fig. llSl with the dashed curve in the upper left panel. These may be divided 
by 2.6 to return to the original local value. J05 actually used a constant- 
extinction correction to Ha to be more consistent with their higher redshift 
correction to [Oil], but the fact that the [Oil] and Ha SFRDs are such 
intrinsically different shapes is by far the dominant effect. 



in the range 8.5 <log(M*/M )< 9.5 down to a limiting SFR of 
~ 0.3M Q yr' 1 . 

We examine the [O II] luminosity function and find a faint 
end slope of af aint ~ —1.5, comparable to the local measurement 
from SDSS. This matches on well to the more [O II]-luminous 
(higher mass) sample from the Gemini Deep Deep Survey (GDDS) 
and ESO public spectroscopy in the GOODS-S field. 

By carefully combining the ROLES data with other spectro- 
scopic surveys, we study the mass dependence of the star forma- 
tion rate density (SFRD) at z~l using two independent SFR in- 
dicators: empirically-corrected [O II] (using the G10 prescription), 
and SED-fitting, which is largely based on rest-frame UV luminos- 
ity (for galaxies with spectroscopic redshifts). By comparing with 
the local SFRD from Stripe 82 of the SDSS (G10), we find that 
both indicators show that the SFRD decreases equally for galaxies 
of all masses (8.5<log(M*/M )<9.5) between z~l and z~0.1. 
The exact change in normalisation depends on the indicator used, 
with the [O II]-based estimate showing a change of a factor of «2.6 
and the UV-based a factor of ~6. We discuss possible reasons for 
the discrepancy in normalisation between the indicators, but note 
that the magnitude of this uncertainty is comparable to the discrep- 
ancy between indicators seen in other z~l works. Our result that 
the shape of the SFRD as a function of stellar mass (and hence 
the mass range of galaxies dominating the SFRD) does not evolve 
between z~l and z~0.1 is robust to the choice of indicator. 

The term cosmic downsizing has been employed to describe 
various different aspects of the cosmic star-formation history as a 
function of galactic mass. In perhaps the most common usage, one 
might expect that the downsizing picture would describe a situation 
in which the peak of the SFRD as a function of mass shifts towards 
lower masses at lower redshifts. This is contrary to the scenario 
seen here, in which the shape of the SFRD does not evolve and 
the normalisation changes equally for all galaxy masses. However, 
it is worth keeping in mind that, considering galax ies with high 
SFRs , these are primarily high mass galaxies (e.g., [Noeske et al.l 
120071) and these will indeed both have dominated the SFRD and 
been more numerous at z~l. This could be considered 'downsiz- 
ing'. In order to avoid the ambiguities in the term 'downsizing', 
we prefer to simply talk in terms of the shape of the SFRD-mass 
here. In the next paper in the series, we will examine the SFRs and 
stellar masses of individual galaxies and confront our observations 
with the latest theoretical models of galaxy formation. 



ACKNOWLEDGMENTS 

This paper includes data gathered with the 6.5-m Magellan Tele- 
scopes located at Las Campanas Observatory, Chile. We thank LCO 
and the OCIW for the allocation of time to this project as part 
of the LDSS3 instrument project. We thank T. Dahlen, H. Flo- 
res, S. luneau, B. Mobasher, C. Ravikumar, E. Vanzella and S. 
Wuyts for graciously providing data and for useful discussions, and 
Hsiao- Wen Chen and Sandra Savaglio for helpful comments. Karl 
Glazebrook and I-H Li acknowledge financial support from Aus- 
tralian Research Council (ARC) Discovery Project DP0774469. 
Karl Glazebrook and Ivan Baldry acknowledge support from the 
David and Lucille Packard Foundation. MLB acknowledges sup- 
port from the province of Ontario in the form of an Early Re- 
searcher Award. 



20 D. G. Gilbank et al. 



REFERENCES 985 

», , , „™„ .T Zhu, G., Moustakas, J., & Blanton, M. R. 2008, ApJ, 701, 86 

Abraham, R. G. et al. 2004, AJ, 127, 2455 

Baldry, I. K. & Glazebrook, K. 2003, ApJ, 593, 258 

Bouwens, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2007, 

ApJ, 670, 928 

Brinchmann, J., Chariot, S., White, S. D. M., Tremonti, C, Kauff- 
mann, G, Heckman, T., & Brinkmann, J. 2004, MNRAS, 351, 
1151 

Calzetti, D., Armus, L., Bohlin, R. C, Kinney, A. L., Koornneef, 

J., & Storchi-Bergmann, T. 2000, ApJ, 533, 682 
Conselice, C. J. et al. 2007, MNRAS, 381, 962 
Cowie, L. L. & Barger, A. J. 2008, ApJ, 686, 72 
Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 

112, 839 

Damen, M., Labbe, I., Franx, M., van Dokkum, P. G, Taylor, 

E. N., & Gawiser, E. J. 2009, ApJ, 690, 937 
Dave, R. 2008, MNRAS, 385, 147 
Davies, G. T. et al. 2009, MNRAS, 395, L76 
Dunne, L. et al. 2009, MNRAS, 394, 3 
Fioc, M. and Rocca-Volmerange, B. 1997, A&A, 326, 950 
Forster Schreiber, N. M. et al. 2006, AJ, 131, 1891 
Gehrels, N. 1986, ApJ, 303, 336 

Gilbank, D. G, Baldry, I. K., Balogh, M. L., Glazebrook, K., & 

Bower, R. G. 2010, MNRAS, submitted 
Glazebrook, K. et al. 2004a, Nature, 430, 181 
Glazebrook, K. & Bland-Hawthorn, J. 2001, PASP, 113, 197 
Glazebrook, K., Peacock, J. A., Miller, L., & Collins, C. A. 1995, 

MNRAS, 275, 169 
Glazebrook, K., Tober, J., Thomson, S., Bland-Hawthorn, J., & 

Abraham, R. 2004b, AJ, 128, 2652 
Hogg, D. W. 1999, ArXiv Astrophysics e-prints 
Hopkins, A. M. & Beacom, J. F. 2006, ApJ, 651, 142 
Juneau, S. et al. 2005, ApJ, 619, L135 
Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 
Kewley, L. J., Geller, M. J., & Jansen, R. A. 2004, AJ, 127, 2002 
Lacey, C. G, Baugh, C. M., Frenk, C. S., Silva, L., Granato, G. L., 

& Bressan, A. 2008, MNRAS, 385, 1155 
Lilly, S. J., Le Fevre, O., Hammer, F, & Crampton, D. 1996, ApJ, 

460, L1+ 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., 

Steidel, C. C, & Fruchter, A. 1996, MNRAS, 283, 1388 
Mobasher, B. et al. 2009, ApJ, 690, 1074 

Moustakas, J., Kennicutt, Jr., R. C, & Tremonti, C. A. 2006, ApJ, 
642, 775 

Noeske, K. G. et al. 2007, ApJ, 660, L43 
Pannella, M. et al. 2009, ApJ, 698, LI 16 
Reddy, N. A. & Steidel, C. C. 2009, ApJ, 692, 778 
Salpeter, E. E. 1955, ApJ, 121, 161 
Savaglio, S. et al. 2005, ApJ, 635, 260 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 
525 

Somerville, R. S., Lee, K., Ferguson, H. C, Gardner, J. P., Mous- 
takas, L. A., & Giavalisco, M. 2004, ApJ, 600, L171 

Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, M., 
& Pettini, M. 1999, ApJ, 519, 1 

Tran, K.-V. H., Franx, M., Illingworth, G. D., van Dokkum, P., 
Kelson, D. D., Blakeslee, J. P., & Postman, M. 2007, ApJ, 661, 
750 

Vanzella, E., et al. 2008, A&A, 478, 83 

Wuyts, S., Labbe, I., Forster Schreiber, N. M., Franx, M., Rud- 
nick, G, Brammer, G. B., & van Dokkum, P. G. 2008, ApJ, 682, 



